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Abstract. The discrepancy between the oxygen abundances in high-metallicity Hu regions determined through 
the T e - method (and/or through the corresponding "strong lines - oxygen abundance” calibration) and that 
determined through the model fitting (and/or through the corresponding "strong lines - oxygen abundance” 
calibration) is discussed. It is suggested to use the interstellar oxygen abundance in the solar vicinity, derived 
with very high precision from the high-resolution observations of the weak interstellar OIA1356 absorption lines 
towards the stars, as a "Rosetta stone” to verify the validity of the oxygen abundances derived in Hu regions 
with the T e - method at high abundances. The agreement between the value of the oxygen abundance at the 
solar galactocentric distance traced by the abundances derived in H II regions through the T e - method and that 
derived from the interstellar absorption lines towards the stars is strong evidence in favor of that i) the two-zone 
model for T e seems to be a realistic interpretation of the temperature structure within Hu regions, and ii) the 
classic T e - method provides accurate oxygen abundances in Hu regions. It has been concluded that the "strong 
lines - oxygen abundance” calibrations must be based on the Hu regions with the oxygen abundances derived 
with the T e - method but not on the existing grids of the models for H II regions. 
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1. Introduction 

An investigation of variations of chemical properties 
among galaxies is very important for the development 
of the theory of the structure and evolution of galax¬ 
ies. Accurate abundances are necessary for such investi¬ 
gations. Good spectrophotometry of Hu regions is now 
available for a large number of galaxies, and the realibil- 
ity of abundances is mainly defined by the method for 
abundance determination in Hu regions. 

Abundance in Hu regions can be derived from 
measurements of temperature-sensitive line ratios, such 
as [OIII]AA4959,5007/[OIII]A4363. Following Stasinska 
(2002b), this classical T e - method will be referred to as 
the direct empirical method. The abundance in HII re¬ 
gions can be also derived through photoionization model 
fitting. This method for abundance determination will be 
referred to as the theoretical (or model) method. 

Unfortunately, in oxygen-rich Hil regions the 
temperature-sensitive lines such as [OIIIJA4363 are often 
too weak to be detected. For such H n regions, abundance 
indicators based on more readily observable lines were 
suggested (Pagel et al. 1979; Alloin et al. 1979). The 
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oxygen abundance indicator R 2 3 = ([Oil]AA3727,3729 
+ [OIII]AA4959,5007)/H / 3, suggested by Pagel et al. 
(1979), has found widespread acceptance and use for the 
oxygen abundance determination in HII regions where 
the temperature-sensitive lines are undetectable. The 
strategy of this way of abundance determination is very 
simple: the relation between strong oxygen line intensities 
and oxygen abundances is established based on the Hil 
regions in which the oxygen abundances are determined 
through the T e - method, and then this relation is 
used for the abundance determination in Hu regions in 
which the temperature-sensitive lines are not available. 
The relation (between strong oxygen line intensities 
and oxygen abundances) established on the basis of HII 
regions in which the oxygen abundances are determined 
through the T e - method (direct empirical method) will 
be referred to as empirical calibration. 

The grids of photoionization models are often used to 
establish the relation between strong oxygen line inten¬ 
sities and oxygen abundances (Edmunds & Pagel 1984, 
McCall et al. 1985, Dopita & Evans 1986, Kobulnicky et 
al 1999, Kewley & Dopita 2002, among others). The re¬ 
lation (between strong oxygen line intensities and oxygen 
abundances) established on the basis of the grids of the 
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photoionization models for HII regions will be referred to 
as theoretical or model calibration. 


The early calibrations were one-dimensional (Edmunds 
& Pagel 1984, McCall et al. 1985, Dopita & Evans 1986, 
Zaritsky et al. 1994), i.e. the relation of the type O/H = 
/(i? 23 ) was used. It has been shown (Pilyugin 2000, 
2001 a,b) that the error in the oxygen abundance derived 
with the one-dimensional calibrations involves a system¬ 
atic error. The origin of this systematic error is evident. 
In a general case, the intensities of oxygen emission lines 
in spectra of HII region depend not only on the oxygen 
abundance but also on the physical conditions (hardness 
of the ionizing radiation and geometrical factors). Then 
in determininig the oxygen abundance from line intensi¬ 
ties the physical conditions in the Hii region should be 
taken into account. In the T e - method this is done via 
T e . In one-dimensional calibrations the physical condi¬ 
tions in Hii regions are ignored. Starting from the idea 
of McGaugh (1991) that the strong oxygen lines contain 
the necessary information to determine accurate abun¬ 
dances in (low-metallicity) Hu regions, it has been shown 
(Pilyugin 2000, 2001a,b) that the physical conditions in 
HII regions can be estimated and taken into account via 
the excitation parameter P. A two-dimensional or para¬ 
metric calibration (P method) has been suggested. A 
more general relation of the type O/H = f(P,R 2 3 ) is 
used in the P - method, compared to the relation of the 
type O/H = /(I? 23 ) used in one-dimensional calibrations. 


It should be stressed that ’’strong lines - oxygen 
abundance” calibrations do not form a uniform family. 
One should clearly recognize that there are two different 
types of calibrations. The calibrations of the first type 
are the empirical calibrations (established on the basis 
of Hii regions in which the oxygen abundances are de¬ 
termined through the T e - method). Two-dimensional 
empirical calibrations both at low and at high metallic- 
ities were recently derived by Pilyugin (2000, 2001a,c). 
The calibrations of the second type are the theoretical (or 
model) calibrations (established on the basis of the grids 
of the photoionization models of Hii regions). The two- 
dimensional theoretical calibrations were recently sug¬ 
gested by Kobulnicky et al. (1999) and Kewley & Dopita 
( 2002 ). 


Thus, at the present day there actually exist two scales 
of oxygen abundances in Hu regions. The first (empirical) 
scale corresponds to the oxygen abundances derived with 
the T e method or with empirical calibrations (the P 
method). The second (theoretical or model) scale cor¬ 
responds to the oxygen abundances derived through the 
model fitting or with theoretical (model) calibrations. The 
comparison of these scales of oxygen abundances in Hu 
regions and their evaluation are the goals of the present 
study. 


2. Oxygen abundances in H 11 regions: model 
fitting versus T e - method 

Fig.[l] shows the two-dimensional empirical calibration ob¬ 
tained by Pilyugin (2000, 2001c) (low-metallicity range) 
and Pilyugin (2001a) (high-metallicity range). Each O/H 
R 23 relation is labeled with corresponding value of the 
excitation parameter P which is defined as 

[0/77] AA4959,5007 

" [Oil] A3727 + [07/7]AA4959,5007' ‘ ' ’ 

The points in Fig.[l] are Hu regions with oxygen abun¬ 
dances determined through the T e - method (compilation 
of data from Pilyugin 2000, 2001a). 

Fig.|| shows the two-dimensional theoretical calibra¬ 
tion reported by Kobulnicky et al. (1999). This calibra¬ 
tion is based on the grid of the models of Hii regions 
after McGaugh (1991). The O/H = /(R 23 ) relations af¬ 
ter Kobulnicky et al. (1999) are shown in Fig.|| with lines. 
Every line is labeled with corresponding value of the exci¬ 
tation parameter P. The ionization parameter y was used 
in the two-dimensional theoretical calibration reported by 
Kobulnicky et al. (1999). The ionization parameter y was 
defined as 


, [0777JAA4959,5007 

9 = 108 -[0//]A3727 


( 2 ) 


As can be seen from Eq. ( 0 ) and Eq. (|^), the excitation pa¬ 
rameter P used in the two-dimensional empirical calibra¬ 
tion of Pilyugin (2000, 2001a,c) and the ionization param¬ 
eter y used in the two-dimensional theoretical calibration 
reported by Kobulnicky et al. (1999) are connected by a 
simple relation 


V = log 


P 

1 -P’ 


( 3 ) 


McGaugh (1991) has concluded that his grid of H 11 re¬ 
gion models agrees with existing oxygen abundances de¬ 
termined through the T e - method at low metallicities 
(his Figure 14), the rms of the residuals between his abun¬ 
dances and published abundances based on temperature- 
sensitive line ratios is 0.05dex with no zero-point offset. 
Inspection of Fig.|] confirms that the theoretical calibra¬ 
tion of Kobulnicky et al. (1999), based on McGaugh’s grid 
of H 11 region models, agrees quantitatively with the re¬ 
cent oxygen abundances derived through the direct em¬ 
pirical method (the T e - method) at low metallicities. A 
comparison of Fig.[j] with Fig|| shows that the discrep¬ 
ancy between the theoretical calibration of Kobulnicky 
et al. and the empirical calibration of Pilyugin is neg¬ 
ligible small for the very low-metallicity (12+logO/H 
around 7.3), high-excitation (P around 0.95) Hii regions, 
but the discrepancy increases with increasing metallicity 
and with decreasing excitation parameter, reaching the 
value of AlogO/H around 0.15 dex for Hii regions with 
12+logO/H around 7.9. The agreement between the the¬ 
oretical calibration of Kobulnicky et al. and the empirical 
calibration of Pilyugin disappears for H 11 regions that lie 
on the upper branch of the O/H - R 23 diagram. 
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Fig. 1. The two-dimensional empirical calibration ob¬ 
tained by Pilyugin (2000, 2001c) (low-metallicity range) 
and Pilyugin (2001a) (high-metallicity range). Each O/H 
- R 23 relation is labeled with corresponding value of the 
excitation parameter P. The points are H 11 regions with 
oxygen abundances determined through the T e - method 
(compilation of data from Pilyugin 2000, 2001a). 



Fig. 2. The two-dimensional theoretical calibration re¬ 
ported by Kobulnicky et al. (1999). The lines are the O/H 
R 23 relations for different values of the excitation pa¬ 
rameter P. Every line is labeled with corresponding value 
of the excitation parameter P. The points are the same 
data as in Fig.Q. 


It should be noted that the agreement (at least quan¬ 
titative) between the theoretical calibrations and the em¬ 
pirical calibrations at low metallicities is not a general 
rule. For example, there is a significant discrepancy be¬ 
tween two-dimensional theoretical calibration of Kewley & 
Dopita (2002) and the empirical calibration of Pilyugin at 
high and at low metallicities. The theoretical calibration 
of Kewley & Dopita is also in conflict with the theoretical 
calibration of Kobulnicky et al. 

Thus, the theoretical calibration of Kobulnicky et al. 
(1999) agrees quantitatively with the empirical calibra¬ 
tions of Pilyugin (2000, 2001a,b) at low metallicities and 
these calibrations are in conflict at high metallicities. 

3. Verification of the validity of the T e - method 
at high metallicities 

The validity of the T e - method at high metallicities has 
been questioned in a number of investigations. According 
to Stasinska (2002a,b), at high metallicities large tempera¬ 
ture gradients are expected in ionized nebulae. Therefore, 
the T e - method based on [OIIIJA4363/5007 will under¬ 
estimate the abundances of heavy elements, since the 
[OIIIJA4363 line will be essentially emitted in the high 
temperature zones, inducing a strong overestimate of the 
average electron temperature. Therefore, although with 
very large telescopes it will now be possible to measure 
[OIIIJA4363 even in high metallicity giant Hil regions, 
one should refrain from interpreting this line in the usual 
way. Doing this, one would necessary find sub-solar oxygen 
abundances, even for giant H 11 regions with metallicities 
well above solar. 

Thus, the validity of the T e - method has been verified 
by comparison with the Hu region models. As it can be 
seen in the previous section, the recent H 11 region mod¬ 
els are not indisputable even at low metallicities. Why 
should one expect that Hil region models provide more 
realistic abundances compared to the T e - method at high 
metallicities? Indeed, according to Stasinska (2002b) this 
would be true if the constraints were sufficiently numerous 
(not only on emission line ratios, but also on the stellar 
content and on the nebular gas distribution) and if the 
model fit were perfect (with a photoionization code treat¬ 
ing correctly all the relevant physical processes and using 
accurate atomic data). These conditions are never met in 
practice. Abundances are not necessary better determined 
from model fitting. 

Then, the validity of the T e - method at high metallici¬ 
ties cannot be indisputably confirmed or rejected by com¬ 
parison with the recent HII region models. Fortunately, 
there is another way to verify the validity of the T e - 
method at high metallicities. 

High-resolution observations of the weak interstellar 
OIA1356 absorption lines towards the stars allow one to 
determine interstellar oxygen abundance in the solar vicin¬ 
ity with very high precision. It should be noted that this 
method is in fact model-independent. These observations 
yield a mean interstellar oxygen abundance of 319 O atoms 
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per 10 6 H atoms (or 12+log(0/H) = 8.50) (Meyer et al. 
1998, Sofia & Meyer 2001). There are no statistically sig¬ 
nificant variations in the measured oxygen abundances 
from line of sight to line of sight; the rms scatter value 
for these oxygen abundances is low, ±0.05dex. Out to 1.5 
kpc, the oxygen abundances are stable in diffuse clouds 
with different physical conditions as measured by the frac¬ 
tion of H in the form of H 2 . 

Caplan et al (2000) and Deharveng et al (2000) have 
analysed Galactic Hii regions and have obtained the 
slope -0.0395 dex/kpc with central oxygen abundance 
12+log(0/H) = 8.82 and 12+log(0/H) = 8.48 at the so¬ 
lar galactocentric distance. All the available spectra of 
Galactic H 11 regions with measured [OIIIJA4363 lines were 
compiled by Pilyugin et al (2002b), and oxygen abun¬ 
dances in Galactic Hu regions were recomputed in the 
same way, using the T e - method. These data result 
in oxygen abundance 12+log(0/H) = 8.50 at the solar 
galactocentric distance although the dispersion in derived 
abundances is relatively large. Thus, the value of the oxy¬ 
gen abundance at the solar galactocentric distance derived 
from consideration of the H 11 regions is in agreement with 
that derived with high precision from the interstellar ab¬ 
sorption lines towards the stars. The agreement between 
the value of the oxygen abundance at the solar galactocen¬ 
tric distance traced by the abundances derived in H 11 re¬ 
gions through the T e - method and that derived from the 
interstellar absorption lines towards the stars is strong ev¬ 
idence in favor of that i) the two-zone model for T e seems 
to be a realistic interpretation of the temperature struc¬ 
ture within Hii regions, and ii) the classic T e - method 
provides accurate oxygen abundances in H 11 regions up to 
oxygen abundances as large as 12+log(0/H) = 8.604-8.70. 
Thus, one can concluded that the Hii regions with T e - 
abundances provide a more reliable basis for calibration 
than the Hu region models. 

Unfortunately, there are no H 11 regions with mea¬ 
sured [OIIIJA4363 lines at highest metallicities (or 
at lowest logR. 23 , Fig.|l]|). There are only two high- 
metallicity, 12+log(0/H) ~ 8.9, Hu regions with mea¬ 
sured temperature-sensitive lines, see Kinkel & Rosa 
(1994), Castellanos et al. (2002). Therefore, we have to 
use an extrapolation of the calibration to derive the abun¬ 
dances in Hii regions with lowest value of logR 23 - The 
oxygen abundances in those H 11 regions predicted by the 
empirical calibration are significantly lower compared to 
abundances predicted by theoretical (model) calibrations. 
There is hovewer indirect way to test the reality of these 
predictions. The present-day oxygen abundance in the so¬ 
lar vicinity is 12+log(O/H)=8.50, and for the present-day 
gas mass fraction, 0.154-0.20 appears to be a reasonable 
value (Malinie et al. 1993). Fig.|| shows the O/H - /r dia¬ 
gram, where /i is the gas mass fraction. The prediction of 
the closed-box model for the chemical evolution of galax¬ 
ies is presented by the solid line in Fig.|j. The constant 
for O/H values was chosen in such manner that the value 
of the gas mass fraction n=0.2 corresponds to the oxy¬ 
gen abundance as large as 12+logO/H=8.5 as in the so- 



Fig. 3. The oxygen abundance versus gas mass fraction 
diagram. The constant was chosen in such a way that 
the value of the gas mass fraction fi=0.2 corresponds to 
the oxygen abundance as large as 12+logO/H=8.5. This 
comes from the data for the solar vicinity interstellar 
medium (see text). 


lar vicinity (square in Figj^). Inspection of Fig.^j shows 
that there is no the possibility of large central (inter¬ 
sect) oxygen abundances (as large as 12+logO/H ~ 9.50 
and higher) derived with theoretical (model) calibrations 
(Zaritsky et al 1994, Garnett et al 1997). On the contrary, 
the low central oxygen abundances predicted by the em¬ 
pirical calibration (Pilyugin et al. 2002a) fit this picture 
well. 

It has been known for a long time that permitted lines 
in HII regions indicate higher oxygen abundances than for¬ 
bidden lines. Peimbert (1967) suggested that the presence 
of spatial temperature fluctuations in gaseous nebulae can 
significantly influence the oxygen abundance in Hii re¬ 
gions derived from forbidden lines. In contrast, permitted 
lines are almost independent of such variations and, in 
principle, they should be more precise indicators of the 
true chemical abundances. Esteban et al. (1998) found 
oxygen abundances for two positions in the Orion neb¬ 
ula from permitted and forbidden lines. The oxygen abun¬ 
dances derived from forbidden lines are coincident for both 
positions in the Orion nebula (12 + logO/H = 8.47) and 
agree well with the interstellar oxygen abundance in the 
solar vicinity derived from interstellar absorption lines to¬ 
wards the stars (12 + logO/H = 8.50). They found from 
permitted lines the oxygen abundance 12 + logO/H = 
8.61 for position 1 and 12 + logO/H = 8.68 for position 
2 , although the abundances obtained from the different 
multiplets observed show significant dispersion. The gas- 
phase oxygen abundance 12 + logO/H = 8.64 and the to¬ 
tal (gas+dust) oxygen abundance as large as 12+logO/H 
= 8.72 were proposed by Esteban et al. (1998) for the 
Orion nebula. If permitted lines indicate the true oxy¬ 
gen abundance in the Orion nebula, then the uncertainty 
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around 0.1 dex in the oxygen abundances derived from 
forbidden lines cannot be excluded. However, the origin of 
the discrepancy between abundances derived from permit¬ 
ted and forbidden lines is not indisputable (see discussion 
in Stasinska, 2002b). If the total (gas+dust) oxygen abun¬ 
dance in the Orion nebula coincides with the total oxygen 
abundance in the interstellar medium in the solar vicinity, 
and if permitted lines provide true gas-phase oxygen abun¬ 
dance in the Orion nebula, then the difference between 
gas-phase oxygen abundance in the Orion nebula and in 
the interstellar gas suggests that the depletion of oxygen 
into dust grains in the interstellar medium is around 0.1 
dex higher than in the Orion nebula (and absolute deple¬ 
tion of oxygen in the interstellar medium is around 0.2 
dex). It can be verified, in principle, by measurement of 
the abundance of the same noble gas in the interstellar 
medium and in the Orion nebula. Unfortunately, measure¬ 
ments of the abundance of only the noble gas krypton in 
the interstellar medium are available (Cardelli & Meyer 
1997), but there is no data on the krypton abundance in 
the Orion nebula. 

Thus it appears that the classic T e - method provides 
more accurate oxygen abundances in H II regions at high 
metallicities as compared to the model fitting, although 
the uncertainty around O.ldex cannot be excluded. As a 
consequence, the empirical calibration appears to be more 
justified than the theoretical (model) calibration. 


4. Conclusions 

The oxygen abundances in Hil regions determined 
through the direct empirical method (the classic T e 
- method) and/or through the corresponding empirical 
’’strong lines - oxygen abundance” calibration are com¬ 
pared with abundances determined through the model fit¬ 
ting and/or through the corresponding theoretical (model) 
’’strong lines - oxygen abundance” calibration. It was 
shown that the theoretical calibration of Kobulnicky et 
al. (1999) agrees quantitatively with the empirical cali¬ 
brations of Pilyugin (2000, 2001a,b) at low metallicities 
and these calibrations are in conflict at high metallicities. 

It is suggested to use the interstellar oxygen abundance 
in the solar vicinity, derived with very high precision from 
the high-resolution observations of the weak interstellar 
OIA1356 absorption line towards the stars, as a ’’Rosetta 
stone” to verify the validity of the oxygen abundances de¬ 
rived in high-metallicity H II regions with the T e - method. 
The agreement between the value of the oxygen abundance 
at the solar galactocentric distance traced by the abun¬ 
dances derived in Hil regions through the T e - method 
and that derived from the interstellar absorption lines to¬ 
wards the stars is strong evidence in favor of that i) the 
two-zone model for T e seems to be a realistic interpre¬ 
tation of the temperature structure within Hu regions, 
and ii) the classic T e - method provides accurate oxy¬ 
gen abundances in Hil regions, although the uncertainty 
around O.ldex cannot be excluded. 


It has been concluded that at high metallicities the 
’’strong lines - oxygen abundance” calibrations must be 
based on the HII regions with the oxygen abundances de¬ 
rived through the T e - method but not on the existing 
grids of the models for Hu regions. 
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